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We show that long-ranged superconducting order is not necessary to guarantee the existence of Majorana 
fermion zero modes at the ends of a quantum wire. We formulate a concrete model which applies, for instance, 
to a semiconducting quantum wire with strong spin-orbit coupling and Zeeman splitting coupled to a wire with 
algebraically-decaying superconducting fluctuations. We solve this model by bosonization and show that it 
supports Majorana fermion zero modes. We show that electron backscattering in the superconducting wire, 
which is caused by potential variations at the Fermi wavevector, generates quantum phase slips which cause 
a splitting of the topological degeneracy which decays as a power law of the length of the superconducting 
wire. The power is proportional to the number of channels in the superconducting wire. Other perturbations 
give contributions to the splitting which decay exponentially with the length of either the superconducting or 
semiconducting wires. We argue that our results are generic and apply to a large class of models. We discuss 
the implications for experiments on spin-orbit coupled nanowires coated with superconducting film and for 
LaA10s/SrTi03 interfaces. 

PACS numbers: 



I. INTRODUCTION 

Kitaev 1 showed that a class of superconducting quantum 
wires supports a pair of Majorana fermion zero modes, one 
at each end. Lutchyn et and Oreg et al? discovered that, 
in the presence of a parallel magnetic field, semiconducting 
wires with strong spin-orbit coupling fall in this class if su- 
perconductivity is induced by proximity to a bulk 3D super- 
conductor, see Fig.[T] As a result of the Majorana zero modes, 
the ground state is doubly degenerate. The two states differ 
by fermion parity, which is not locally measurable; therefore, 
they form a protected qubit. Networks of such semiconduct- 
ing wires have been proposed for topological quantum infor- 
mation processing^. 

Long-ranged superconducting order is an essential feature 
of these analyses. Though such order is sufficient, it does not 
seem necessary. Protected Majorana zero modes also exist in 
models of the 5/2 fractional quantum Hall state 8 13 and in Ki- 
taev's honeycomb lattice spin model 14 , and neither of these 
systems has long-ranged superconducting order. Therefore, 
one might expect that a quantum wire with strong supercon- 
ducting fluctuations but no long-ranged order could also sup- 
port Majorana fermion zero modes. Consider, on the other 
hand, a spinless one-dimensional Luttinger liquid, which has 
algebraic order, i.e. the two-point correlation function of the 
superconducting order parameter decays to zero as a power 
of the separation, rather than approaching a constant. Such 
a system has gapless bulk fermionic excitations, so if Majo- 
rana fermion zero modes were found at the ends of such a 
model, there would be nothing protecting them against small 
perturbations. Furthermore, in the absence of superconduct- 
ing order, the two states of a pair of Majorana fermion zero 
modes would have different electric charges - and not merely 
fermion parity. Simply changing the electrostatic potential 
should cause an energy splitting between states with differ- 
ent electric charges. Therefore, one might, instead, conclude 
that long-ranged superconducting order is necessary to protect 



Majorana fermion zero modes in quantum wires. 

In this paper, we show that this is not the case. We con- 
struct a model of a spin-orbit coupled semiconducting wire 
in a magnetic field which is coupled to an s-wave supercon- 
ducting wire with power-law order. The schematic plot of 
the heterostructure is depicted in Fig. [2] We show that this 
model supports Majorana fermion zero modes at the ends of 
the wire. However, a single wire does not support a qubit; at 
least two wires are needed. The basic idea is simple. Con- 
sider Kitaev model of a superconducting quantum wire of 
spinless fermions. 

H = -t Yj( c Ui c i + 4 c i+i) 
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Here, (j) is the phase of the superconducting order parameter. 
Let us assume, for the moment, that is a constant, as in 
Kitaev 's original papei^ and in Refs. 12131 If we rotate the 
fermion operators to the local value of the phase of the order 
parameter: q — ^ e^/ 2 Q, then the Hamiltonian takes the form 

H = -^(cl +1 c,+cJc, +1 ) + |A| ^(qq +1 -cJcJ +1 ) (2) 

i i 

At the special point t = |A|, this Hamiltonian can be di- 
agonalized by introducing the Majorana fermion operators 

72i-i = c» + cj, 72i = (c* - c\)/i\ 

tf = z|A|^ 72z72z+1 (3) 

i 

These operators satisfy ji = 7 J and {7^,7^} = 25ij. Note 
that 71 and j2N do not appear in the Hamiltonian. Therefore, 
the ground state is doubly degenerate: 271727V can be either 
±1 while ^72i72i+i = —1 for 1 < i < N — 1. Apart from the 
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degeneracy of the ground state, there is a gap 2 1 A | to excita- 
tions. The operators 71 and 72 at are Majorana fermion zero 
modes, and the qubit which they form, 271727V = =1=1, is pro- 
tected since the two states are distinguished only by fermion 
parity, which cannot be measured by a local operation. Only 
an operator which acts on both sites 1 and N can affect it. 
Away from the special point, t = |A|, the physics is very 
similar: there is a gap in the bulk above two nearly degener- 
ate ground states which have an energy splitting ~ e -Na/£^ 
where £ is inversely proportional to the bulk gap and a is 
the lattice spacing. This phase persists to the more physi- 
cal |A| <C t limit. Electron-electron interactions in the wire 
determine the region of the phase diagram occupied by this 
phase™. 

Now suppose that <p is a fluctuating dynamical field. We can 
still perform a change of variables similar to the one which we 
made in going from Eq. ([]} to Eq. ([2]). This will remove the 
phase of the order parameter from the second term in Eq. ([T]), 
the pairing term. However, it will introduce a coupling be- 
tween the fermions and gradients of the order parameter. If 
these terms can be neglected, then we will have mapped a 
model with fluctuating order parameter to one with fixed or- 
der parameter which is decoupled from the fluctuations of <p\ 
therefore, it will have Majorana fermion zero modes. How- 
ever, there are some subtleties involved in the change of vari- 
ables from Ci to at when <\> fluctuates. These are most easily 
handled using a bosonized formulation of the electronic de- 
grees of freedom in the wire. We find a special point in Sec- 
tions III and IV at which the bosonized formulation simplifies 



and allows us to completely analyze the model. We then show 
in Section [V] that our analysis is qualitatively unchanged by 
perturbations which take the system away from the special 
point. 

The technical subtleties alluded to above have a physical 
origin related to the conservation of charge. Note that the 
ground state energy has the form 



E(N) = NS + E i 
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for even and odd electron numbers N, respectively. (See Ref. 
[19] for the analogous relation for paired quantum Hall states.) 
The signature of Majorana fermion zero modes at the end- 
points of a wire is that i^odd = E even . In a superconducting 
system without zero modes, we would have E Q ^ > E even . 
The difference E Q ^ — ^even would simply be the energy cost 
of an unpaired electron. In the presence of zero modes, this 
cost vanishes. As may be seen from ([4]), however, a single 
wire does not have degenerate states unless the electrostatic 
potential is tuned so that £ = 0. 

If, however, we consider two such wires, then there are 
two degenerate states for fixed total electron number with- 
out any fine-tuning. Suppose that there are 2N electrons in 
the system. Let us denote the energy of the two wires, iso- 
lated from each other, by E\(N), E^iV). They are given by 
lil with £ ( 2 ) and 25<J> 25® odd taking the place 



4dd = ^evL Then 

Ei(N) + E 2 (N) = Ei(N -m) + E 2 (N + m) (5) 

for any m, so long as £^ = £^ 2 \ Now suppose that the 
two semiconducting wires are coupled to the same (power- 
law) 5-wave superconducting wire (which is assumed to be 
much longer than either semiconducting wire so that it can 
be coupled to both while keeping them far apart), so that the 
electrochemical potential must be the same in the two wires. 
Then £^ = £ ^ . Furthermore, Cooper pairs can tunnel from 
either semiconducting wire to the superconductor. Therefore, 
rather than a degenerate ground state for each value of m in 
([5]), there will be two nearly degenerate states, corresponding 
to an even or odd number of electrons in each wire. Such a 
protected qubit exists for any fixed electron number. If the 
electron number were odd, then the two states would corre- 
spond, instead, to (a) even electron number in wire 1 , odd in 
wire 2; and (b) odd electron number in wire 1, even in wire 2. 

These arguments are supported by explicit calculations in 
Sections [TTTJ and |lVj First, we show in Section [TT] how the 
topological degeneracy is manifested when a semiconducting 
nanowire is coupled to a bulk 3D superconductor. Pair tun- 
neling between the wire and the 3D superconductor is repre- 
sented by a term in the bosonized effective Hamiltonian of the 
form 



tun . oc sin 29 



(6) 



where 6 is the bosonic field satisfying p = ^d t 0, where p is 
the charge density. The two ground states of the system cor- 
respond to the two minima of sin 20 as a function of 6. As we 
discuss in Section [TTJ these two states differ in fermion par- 
ity, as expected for a pair of Majorana zero modes. Further- 
more, if the two ends of the wire are connected to form a ring, 
then the ground state degeneracy disappears because only the 
equal amplitude superposition of the two minima is allowed 
for periodic boundary conditions of the electrons (while the 
orthogonal superposition occurs for anti-periodic electronic 
boundary conditions). When we turn in Sections|lTl|and|Tv|to 



of £ and i£even,odd- If there are Majorana zero modes at the 
endpoints of both wires in isolation, then E^ d = Eill n and 



the case in which the superconductor is also one-dimensional 
and, therefore, does not have long-ranged order, our analysis 
will depend on a careful treatment of the target space of the 
bosonic fields. The periodicity conditions satisfied by these 
fields encode the quantization of charge, and the ground state 
degeneracy cannot be counted properly without accounting 
for them. The use of bosonization techniques also requires 
a careful treatment of locality: putative Majorana modes in a 
transformed system may simply be a reflection of a sponta- 
neously broken global Z 2 in the original variables, c.f. the du- 
ality between the transverse field Ising model and a Majorana 
wire. We wish to stress the topological nature of the Majorana 
degeneracy in our model: no local observable can distinguish 
the two states. A key feature of these models is that there 
is a single-fermion gap even though there are gapless super- 
conducting phase fluctuations, as is already apparent in ^ if 
the second line is benign (as we show it to be). This may be 
viewed as a form of the "spin-gap proximity effect' 2021 . This 
gap protects the Majorana fermion zero modes. However, as 
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we show below, in addition to fermion tunneling events which 
lift the topological degeneracy even in models with long-range 
superconducting order, there is another error-causing process 
involving quantum phase slips which will have a vanishingly 
small probability of occurring in a bulk 3D superconductor. 
The effect of a quantum phase slip in the middle of a super- 
conducting wire can be understood as that of a vortex encir- 
cling a pair of Majorana zero modes. Such a process results in 
reading out the fermionic parity via the Aharonov-Casher ef- 
fect and effectively leads to a splitting of the degeneracy. We 
show that backscattering from impurities generates quantum 
phase slips in the middle of the wire and causes a splitting of 
the topological degeneracy which decays algebraically with 
the size of the system rather than exponentially. However, the 
exponent is proportional to the number of channels in the su- 
perconducting wire. Thus, by making a superconducting wire 
with sufficiently-many channels, we can make the splitting de- 
cay as a high power of the length. 

When this is the case, it is sufficient for the wires and wire 
networks of RefsP^to be in proximity to systems with power- 
law superconducting order; long-ranged order is not neces- 
sary. Consequently, it may be possible to sputter supercon- 
ducting grains onto the semiconducting wire or to coat it with 
superconducting film of a finite thickness. This is important 
because it may be difficult to tune a semiconducting wire be- 
tween topological and non-topological phases by applying a 
gate voltage if it is in contact with a bulk superconductor 
which will presumably fix its chemical potential. 

Recently, it has been shown that quasi- ID wires can be 
"written" on LaAlOa/SrTiOs interfaces 22 which have substan- 
tial Rashba spin-orbit coupling 23 . These wires show strong 
superconducting fluctuations. As we will discuss in detail 
elsewhere 24 , a possible model for this system is a spin-orbit 
coupled quantum wire in contact with local superconduct- 
ing regions which fail to percolate across insulating SrTiC>3 
but can induce still superconducting fluctuations in quantum 
wires at the LAO/STO interface. Our results imply that these 
superconducting fluctuations may be sufficient to support Ma- 
jorana zero modes at the ends of such wires if a parallel mag- 
netic field is applied. 

II. A SEMICONDUCTOR NANOWIRE COUPLED TO A 
BULK 3D SUPERCONDUCTOR: BOSONIZED 
FORMULATION 

Before introducing our model, we briefly review the pro- 
posal for realizing Majorana quantum wires in semiconductor- 
superconductor heterostructures 2 ^ and recast it in bosonic 
form. Its basic ingredient is a semiconductor nanowire with 
strong spin-orbit interactions. Superconductivity is induced 
via the proximity effect. The Hamiltonian for the nanowire is 
(fi = 1): 

rL/2 / q2 \ 

H NW = 1 dx ^t(x)\-—^-fi^iaa y d x +V x a x j 0), 

fL/2 

H P = / dx [A ^t^ + h - c -} • (7) 

J -L/2 




FIG. 1 : A semiconductor nanowire in contact with a bulk 3D super- 
conductor. 

where m*, /i and a are the effective mass, chemical po- 
tential and strength of spin-orbit Rashba interaction, respec- 
tively. An in-plane magnetic field B x leads to a spin splitting 
V x = gsM^BB x /2, where gsM and \ib are the g-factor in the 
semiconductor and the Bohr magneton, respectively. In the 
simplest model for the nanowire, we assume that the semicon- 
ductor nanowire (NW) is in tunneling contact with a bulk 3D 
superconductor (SC), as depicted in Figure [T] Then, electron 
tunneling between the NW and the SC leads to the proximity 
effect described by the Hamiltonian Hp . The superconduct- 
ing pairing potential A is assumed to be a static classical 
field and quantum fluctuations of the superconducting phase 
are neglected. 

The nanowire described by the Hamiltonian Ht = #nw + 
Hp can be driven into a non-trivial topological state by ad- 
justing the chemical potential so that it lies in the gap |/i| < 
\/V x — Aq. Under these conditions the Hamiltonian can be 
projected to the lower band of the two bands which form as 
a result of the combined effect of the spin-orbit coupling and 
magnetic field. The low-energy limit of this Hamiltonian then 
takes the same form as Eq. ^ for low energies E <C ^as- 
suming |A| <Ci 2 . Therefore, the topological superconducting 
phase described by H T harbors Majorana fermion operators 
7l and which are zero modes, up to exponential correc- 
tions, localized about the two endpoints: 

7a = ll , {7a, 7b} = 2£ a6 (8) 

[At, 7a] =0 + O(e- L ^). (9) 

{ 7 l, Vv(*)} ~ e-^ L ^ , { lR A*{x)} ~ e-\ x ~ L ^ 

(10) 

Here, £ is the effective coherence length. The presence of 
these zero-modes leads to topological degeneracy up to an 
exponential splitting energy SE oc e~ L ^. The two nearly- 
degenerate states correspond to the two eigenvalues of 271 72 
and have even and odd fermion-parit} 1 , respectively, which 
can be exploited for topological quantum computation^. 

These results were obtained 1 3 using the properties of the 
free fermion band structure embodied by H T . We now re- 
derive them using a bosonic representation. In later sections, 
we will use this representation to analyze the case when there 
is no long -ranged superconducting order, unlike in Ht. First, 
we bosonize the semiconductor Hamiltonian ([7]). In the helical 
regime corresponding to a large Zeeman gap, i^NW can be 



4 



approximated by projecting the system to the lowest subband 
and writing the field operator Sfr(x) = (#), ^i(x)) as 

*(x)*$-(p F )e ipFX c R (x) + $-(-p F )e- ipFX c L (x) (11) 

where the spinor &_(p F ) = (— e lK ( pF \ l) and 

^(pf) = tan -1 (apF/14). Substituting ([IT]) into #nw, 
the Hamiltonian can be written in terms of the spinless right 
and left-moving fermions cr(x) and cl(x) and eventually 
bosonized using c R / L = 7= — e~ l ^ ±( ^~ e ^'. 



/27va 



H } 



NW 



V 

2^ 



L/2 

dx 

-L/2 
L/2 



ic ] L (x)d x c L (x) - ic ] R (x)d x c R (x) 
/ dx[K{d x 0) 2 + K-\d x <l>) 2 ]. 

J -L/2 



(12) 
(13) 

Here v is the fermion velocity v = Pf{~^ 1 °^ „ ) and 

K the Luttinger parameter for the nanowire. The fields (j) and 
satisfy the canonical commutation relation: 



[d x (j>(x), 0(x')\ = m5(x — x') 



(14) 



The charge density and current near wavevector zero are given 
by p = \d x c\> = \d t Q and j = -\dt<$> = \d x Q. The 
fields (j) and can be interpreted as the phase of the density 
at wavevector 2k f and the pair field, respectively: 



-2i(p(x) 



p2k F (x) = 
^pair(^) = ^(x)^(x) 



2i6(x) 



(15) 



For the Hamiltonian Hnw, in which electron-electron in- 
teractions in the semiconductor have been neglected, K = 1, 
the free-fermion value. However, the bosonic representation 
accommodates short-ranged interactions in the nanowire such 
as 



NW int. 



rL/2 

= u dx^x^^x)^,^)^^^) (16) 

J -L/2 

simply by shifting the value of K and rescaling v. K < 1 for 
repulsive interactions and K > 1 for attractive interactions. 
The bosonic form for Hp in Eq. ([7]) is: 



Hp 



.L/2 
/ 

J-L/2 



dx sin (20) 



(2ira) j_ L/2 
Therefore, can be written in the bosonic form 

H T = [ L/ dx(^- [K(d x 0) 2 + K-\d x ^) 2 } 

J-L/2 \ Z7r 



(17) 



(2na) 



sin (20)1 (18) 



This interaction term, Hp, is relevant unless there are very 
strong repulsive interactions in the nanowire. To be more pre- 
cise, the lowest-order RG equation for the dimensionless cou- 
pling y = 2Apa/v is: 



dy 

dl 



= (2-K~ 1 )y 



(19) 



For non-interacting electrons, K = 1, and even for repul- 
sive interactions up until K = 1/2, this is a relevant per- 
turbation. If y is initially small at short distances, then 
we can use Eq. ([19]) to conclude that y(l) ~ 1 at the 
length scale I = ln(£/ao), where the effective coherence 
length, £, in the semiconducting nanowire is given by £ ~ 
ao(v/2Apao) K /( 2K ~ 1 \ Here, ao is the short-distance cut- 
off, which is the shortest length scale at which the effective 
description ( [T8] ) is valid. We can take it to be the coherence 
length or the Josephson length of the bulk 3D superconductor 
but, at any rate, it must be larger than the Fermi wavelength in 
the semiconducting wire. 

At longer length scales, the field is pinned to the minimum 
of sin(20). Since there are two minima, = — 7r/4, 37r/4, 
there are two degenerate ground states in the L — >> oo limit. 
These two ground states are related to each other by the global 
Z 2 symmetry of the model, —> + 7r. To understand this 
symmetry better, it is helpful to note that the fermion parity 
( — 1)^ can be written in the form 



(-1) 



N F = e i(0(L/2)- 



H-L/2)) 



(20) 



Therefore, using the commutation relation ([14]), we see that 
the fermion parity (— 1)^ generates the symmetry transfor- 
mation — » + 7r. Since the two degenerate ground states 
corresponding to — — 7r/4, 37r/4 are transformed into each 
other by fermion parity, the following quantum superpositions 
are fermion parity eigenstates: 



| even, odd) 



-tt/4> ± |3?r/4) 



(21) 



The ends of the wire are crucial for this qubit. If we were to 
connect the two ends of the wire to form a ring of circumfer- 
ence L, then we would expect only a single ground state, not 
a degenerate pair. To see that this is, indeed the case, consider 
the fermion annihilation operators: 



cr,l(x) 



-i(±<t>-e) 



(22) 



Since p = \d x §, the ring will have even fermion parity if the 
boundary conditions on <p are: 

<j>(x + L) = <j>(x) + 2nir 

for integer n. If the fermions have periodic boundary condi- 
tions, cr,l(x + L) = cr,l(x), then the boundary condition 
on must be 

6(x + L) = 0(x) + 2n7r 

for integer n' . Since constant solutions are allowed for 
this boundary condition on 0, the ground state |even) = 

-^{\— 7r /4) + |3tt/4)^ , which is a linear superposition of 
constant solutions, is allowed in this case. This state has 
even fermion parity ( [20] ), so it is consistent with the bound- 
ary conditions on 0. If the ring has odd fermion parity, how- 
ever, then (j){x + L) = <j>(x) + (2n + 1)tt. Consequently, if 
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the fermions have periodic boundary conditions, the boundary 
condition on must be 0(x + L) = 0(x) + (2n' + l)7r. This 
precludes a constant solution. Therefore, the state |odd) = 
4= ( | — 7r/4) — 1 37r/4) ) , which is odd under fermion parity ( 20 ), 
is not an allowed state if the fermions have periodic boundary 
conditions. As expected, we conclude that there is only a sin- 
gle ground state for a ring, in contrast with a line segment 
which has a doubly degenerate ground state. 

The Majorana fermion zero modes of this system are man- 
ifested on a ring by the presence of a corresponding state 
for anti-periodic boundary conditions on the fermions. If 
cr,l(x + L) = —cr^l{x), then for odd fermion parity, 
(j)(x + L) = <j)(x) + (2n + l)7r, the boundary condition on 
must be 0(x + L) = 0(#) + 2n'7r for integer n'. This bound- 
ary condition allows constant solutions, so the ground state is 
| odd) = -^(|-7r/ 4 ) - l 37r / 4 ))- Therefore, the ground state 
with periodic boundary conditions and the ground state with 
anti-periodic boundary conditions have the same energy den- 
sity and opposite fermion parities. This can already be seen 
in the Kitaev chain. On a line segment, the operators 71 and 
72 at do not appear in the Hamiltonian, as we saw in the intro- 
duction. On a ring with periodic boundary conditions, there is 
a term it^Nji- If the boundary conditions are anti-periodic, 
the term is instead —itj2Nji- The ground state energy is the 
same in both cases, but the ground states differ in fermion par- 
ity, ij2Njl = ±1. 

Returning now to the case of open boundary conditions, we 
observe that, for finite L, these two states are split in energy 
because there are instantons which tunnel between the two 
minima. The Euclidean action in the strong coupling limit is 



S 



=— [ 



dxdr 



(d x ef 



\d r of 



y_ 



sin(20) 



(23) 



The splitting is then given by SE ex Nfe~ s °, where So 
is the action of the Euclidean instanton 6q(x,t) satisfying 
6q(x, —00) = —j,6q(x, 00) = and Nf is a prefactor that 
comes from fluctuations. Clearly the lowest action instanton is 
translationally invariant, at least away from x = — L/2, L/2, 
so the problem reduces to a + 1 dimensional problem, with 
action 



Sqm 



dz 



(d z ef + v{e) 



(24) 



where V{6) = % sin(2#) and z = vt. Following Ref. 



So = - / d0y/2(V(9) - E) 

K J-n/4 



7T C 



where E = —y/^ 2 is the energy of the minimum of the po- 
tential. The splitting then scales like 5E oc exp ^— ^, 
as expected. 

Since 20 changes by 2tt while the phase of the bulk super- 
conductor is unchanged, such an instanton can be interpreted 
roughly as the motion of a vortex between the NW and the 
bulk superconductor. (We say "roughly" because our instan- 
ton is a spatially uniform phase slip, rather than a spatially- 
localized vortex.) Since it causes a transition between the 
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(b) 



FIG. 2: A semiconductor nanowire in contact with a ID supercon- 
ducting wire. The superconducting wire could be a coating which (a) 
completely covers the semiconductor or (b) only covers part of it. 



states I — 7r/4) and |3tt/4), it splits the states |even) and |odd). 
Thus, it can also be interpreted as Majorana fermion tunneling 
between the two ends of the wire. 



III. A SINGLE SEMICONDUCTING NANOWIRE 
COUPLED TO AN ALGEBRAICALLY-ORDERED 
SUPERCONDUCTING WIRE 

We now include the effect of quantum fluctuations by re- 
placing the bulk superconductor in the above proposal with 
an s-wave superconducting wire with power-law order. This 
model preserves the overall U(l) charge symmetry (there is 
no spontaneous U(l) breaking) and allows for the study of 
the topological superconducting phase in the particle number- 
conserving setting. For the sake of concreteness and simplic- 
ity, we will take the Hamiltonian for the superconducting wire 
to be the attractive- U Hubbard model. However, our results 
hold for any spin-gapped system with s-wave superconduct- 
ing fluctuations. 

We use the standard bosonization procedure for spinful 
fermions, with the convention^ that 



-.e ^2 



[(r0 p -0 p )+a(r^ CT )] 



(26) 



where r = ± and a = ± for right/left-moving fermion with 
t / I spin, and a the lattice cutoff. The fields <j) p ^ and P ^ 
satisfy the same commutation relations ( [T4| ). In terms of these 
fields, the Hamiltonian for the superconducting wire can be 
written as 



#sc 



tt{p) 1 rrO) 
r L/2 



(25) ff(*> = V JL [ dx [K p {d x e p f + K-' (d x <P P ) 2 } 

Z7T J_ L /2 



(27) 
(28) 



sc 



2tt 



L/2 
L/2 



L/2 

2\U\ r L/2 



dx [K a (d x 6 a ) 2 + K-\d x ^) 2 } (29) 



(27ra) 



I dx cos(2\/20cr) 

J -L/2 



where vp, a and U are the Fermi velocity, the effective cutoff 
length and the interparticle interaction potential, respectively. 

Tunneling between the superconducting wire and the semi- 
conducting wire can be described using a simple model 
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Hamiltonian 

H t = t 



f L/2 

- J-L/2 



where t is the tunneling amplitude and i/j a and 77^ repre- 
sent fermion annihilation operators in the semiconducting and 
superconducting systems, respectively. Given that single- 
electron tunneling into the superconducting wire is suppressed 
due to the presence of the spin gap E g (see below), the domi- 
nant contribution to the action comes from pair hopping. The 
perturbative expansion in t to second order leads to the fol- 
lowing imaginary-time action 

SpB_ = —t 2 ^^dxdrdxdr (31) 
i/>l(x, r)^_ a (x f , r f )rj a (x, r)rj_ a (x f , r') + h.c. 



2, which includes the case of weakly-attractive interactions in 
the superconducting wire, K p ~ 1, and weakly-repulsive in- 
(30) teractions in the semiconducting wire, K ~ 1. 

This model simplifies significantly at the special point 
vp = v and 2K p = K. At this point, one can diag- 



We now analyze the bosonized action. First, the 
spin field <j) a orders as a result of the last term in 
opening a spin gap E g in the superconducting 
The dual field a is disordered, and its correlation 



Eq. 
wire. 



29 



function decays exponentially (e -h 6a ^ x,T ^e ^ 6><T ^ ' ^ ) ) cr ~ 
a/\/x 2 + (^r) 2 exp[— E g ^r 2 + x 2 /v F }. This allows us to 
simplify the action §T\) and make a local approximation 



^PH < 



(27ra) 



J dr J ' dxsm(V20 p -20 S j (32) 



> Eg 1 . Here the 



valid in the long-time limit |r — r 
Cooper pair hopping amplitude Ap is given by Ap 
t — cxpF = s j m ilarly i me proximity-induced gap in the 

^9 y(apF) 2 +V x 2 

perturbative tunneling limit t <C E g . If the field P were 
pinned (i.e. p = 0), we would recover the model consid- 
ered in Refs. 12131 In the present case, however, overall U(l) 
symmetry is not broken due to the presence of fluctuating field 
P . Henceforth we thus analyze the following effective low- 
energy model 

rL/2 



H M = [ dx [K(d x 9) 2 + K~ 1 (d x cf)) 2 ] 

^ J-L/2 

rL/2 

/ dx[K p (d x p ) 2 + K-\d x <l> p ) 2 } 

J-L/2 



VF_ 

2tt 



(2ira) 



L/2 



dx 



L/2 



sin (V26 p - 20} , (33) 



and study the effect of quantum fluctuations of 6 p on the sta- 
bility of the topological superconducting phase. This model is 
quadratic, except for the interaction Ap. The dimensionless 
coupling y = 2Apa/v has RG equation 



dy 
dl 



(Z-lK^-K^y 



(34) 



For + K~ x > 2, this interaction is irrelevant, and we 

can ignore Cooper pair tunneling between the wires. How- 
ever, inter- wire pair tunneling is relevant for + K~ x < 



onalize the Hamiltonian ( [33] ) by introducing new variables 

(9+ = P /V2 + and 6L = p j\[2 - 0: 



v f L/2 

H=— dx [K p (d x + ) 2 ^K;\dM 2 ] 

Z7T J_ L /2 



2tt 



r L/2 

-L/2 
L/2 

-L/2 



dx [K p (^_) 2 +K p - 1 (^-) 2 ] 



Ap r L/2 



(2ira) 



I 



L/2 



dx sin(26L) . 



(35) 



The first line of this Hamiltonian describes gapless super- 
conducting phase fluctuations. The second and third lines, 
which are decoupled from these gapless fluctuations, are iden- 
tical to the Hamiltonian ( [T8] ) for the proximity effect from 
a bulk 3D superconductor with long-ranged superconducting 
order parameter. At this point, the dimensionless coupling 
y = 2Apa/v has RG equation 



dy 

dl 



(^-K; l )y 



(36) 



Therefore, a fermionic gap Ap ~ -^(Apao/v) 1 ^ 2 K p ) 
opens up as a result of the coupling between the wires. 

The single wire model, however, does not exhibit Majo- 
rana degeneracy without fine tuning of the electrostatic po- 
tential. Semiclassically, this is because the moduli space 
of low energy field configurations (i.e. those where 0_ is 
pinned) has only one connected component. Naively, the 
Ap sm(y/20 p — 20) term might lead one to expect two con- 
nected components, corresponding to the two minima = 
— T* ' — "T • However, these are in fact connected in 
the p , moduli space, see Fig[3] One can interpolate from 
one to the other by winding y/20 p —> \[20 p + 2tt and simul- 
taneously winding half as fast, so that — » + n. ^-re- 
mains pinned throughout the interpolation, but the two vacua 
are exchanged. Therefore, there is no potential barrier; the 
field p /y/2 + is free to fluctuate along a flat direction of 
the potential between these two points. Consequently, there 
is just a single vacuum, not two degenerate states. This re- 
flects the conservation of charge: when p j\[2 + has large 
fluctuations, the total charge is fixed. 

Note that, in the argument above, the two minima were ex- 
changed if we could identify \[20 p = \[20 p + 2tt. Naively, 
these two field values are not equivalent since a shift of V20 p 
by 2tt changes the sign of the fermion according to Eq. [25] 
However, (V20 p , V20 a ) = (y/20 p + 2tt, y/20 a + 2tt). Since 
is fixed, (j is disordered, so there is no energy cost for 
shifting a . Thus, we can treat V20 p as 2tt periodic, rather 
than 47r periodic, and the flat direction of the potential con- 
nects the two putative minima. 
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for such a system takes the form: 



FIG. 3: The moduli space of semiclassical vacua is the torus, which 
is here depicted as a cylinder with the top and bottom edges identi- 
fied. For any fixed P , there are two different semiclassical ground 
state, depcited by the intersection points between the dotted vertical 
line and the line of fixed y/20 p — 20 which winds twice around the 
cylinder. If P is, indeed, fixed, as in Section [il] then there is a tun- 
neling barrier between these two ground states. However, if the total 
charge mode 4= P + can fluctuate, as in Section III then the entire 

line of fixed y/20 p — 20 is the same quantum ground state, and there 
is no degeneracy. 



IV. TWO MAJORANA WIRES 



As discussed in the introduction and as we saw in the pre- 
vious section, if the electron number is fixed then states with 
different electron numbers will not be degenerate without fine- 
tuning. However, if we have two semiconducting wires of 
length I coupled to the same superconducting wire, of length 
L, then there will be two degenerate states of the system for 
any fixed total charge. These states correspond to even or odd 
electron numbers in each semiconducting wire, with a con- 
straint that the sum of the parities of the two wires must equal 
the parity of the total electron number. There need not be 
literally two separate wires. We could instead have a single 
wire similar to the spin-orbit coupled semiconducting wire 
of Eq. [7] In the regions -L/2 < x < -L/2 + £ and 
L/2 — £ < x < L/2, we would need to adjust the chem- 
ical potential so that < \fV 2 
-L/2 + £< x < L/2 



< V£ — Aq and in the region 
' we would need > \/V£ — Aq 
Then, the system would be in a topological (power-law) su- 
perconducting phase for L/2 — £ < \x\ < L/2— and in a non- 
topological phase for —L/2 +£ < x < L/2—£. While we will 
sometimes call the region —L/2-\-£ < x < L/2 — £ the "non- 
topological region", we will usually simply treat the system as 
if there were no wire there the section of non-topological wire 
in this region has a qualitatively similar effect to the absence 
of a wire. 



Let us analyze this setup in more detail. The Hamiltonian 



2 wires 



-L/2+i 



L/2 



2tt 



Api 



sin(V20 p - 20i ) 



(27ra) 

■ / dx[^[K 2 (d x 1 ) 2 ^K^\d x ^ 1 ) 2 } 

JL/2-1 \ Z7r 

sm(V20 p - 20 2) 



A P2 

(27m) 



L/2 



+ / _ dx^ [K p (d x p ) 2 + K-\d x <\> p f 



L/2 



(37) 



The first two lines are the Hamiltonian for the first semicon- 
ducting wire, of length K<I and its Josephson coupling to 
the superconducting wire of length L. The third and fourth 
lines are the analogous terms for the second semiconducting 
wire. The final line reflects the charge degrees of freedom of 
the Hamiltonian for a wire with power-law superconducting 
fluctuations. The gapped spin degrees of freedom have been 
integrated out. In Eq. ( [37] ) we have neglected exponentially 
small corrections ex exp[— E g {L — 21) /vp) due to tunneling 
between the wires, see Sec.|V]for details. 

As in the single wire case, we introduce the fields: 

0+{x) = ^0 p (x)u 12 (x) + e p {x){\ - u 12 (x)) 
+ 9i(x)ui(x) + 2 (x)u 2 (x) 
0-(x) = ^0 p (x)u 12 (x) - O^u^x) - 2 (x)u 2 (x) (38) 

where u\{x) = 1 for —L/2 < x < —L/2 + £ and u\{x) = 
otherwise; u 2 (x) = u(—x); and ui 2 (x) = u\(x) + u 2 (x). 
The field 6-\x) is only defined for L/2- £ < \x\ < L/2. 
Then, for v\ = v 2 — v p = v and K\ = K 2 = 2K p = K, the 
Hamiltonian takes the form: 



2 wires — 
L/2 



'2 / 

j dx(u 12 (x)^ [K p {d x 6-f + K^(d x (t>-) 2 ] 

+ (2^0 ( A f i^i ( x ) + &P2U 2 (x)) sin(26L) 

+ ^ [K p {dJ + ) 2 + K;\d x <P + ) 2 ]^j (39) 

Naively, this Hamiltonian has four semi-classical ground 
states: 0-(x) — <$\U\{x) + (f2U 2 (x), with (pi y2 = 

However, by acting with (—1)^ , (—1)^ , and 

( — 1)^ )+7V f } on any one of these states, we can obtain 
the other three. Thus, we can form two quantum superposi- 
tions of these states with (— 1) n f ) + n f ) = 1 and two with 
(— 1)^ )+Ar F } = —l. If we fix the total electron number, 
then one of these two sets will be allowed. 

The argument which led us to conclude that a single 
nanowire has no ground state degeneracy now shows that the 
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two wire system ( [37] ) at fixed electron number has two (nearly) 
degenerate ground states. There are two connected compo- 
nents in the moduli space of low energy field configurations. 
Here there is no electrostatic potential breaking the expo- 
nentially small degeneracy; the leading contributions instead 
come from instantons, as we shall see below. To see this, note 



first that 0i and 02 can be pinned to either -g= — ^ or -g= — ^ , 
naively leading to four semiclassical ground states. However, 
as above, one can wind \[20 p — » \[20 p + 2tt, thus connect- 
ing the ground state (tp u <p 2 ) = (j,j) with (^f , ), and 
( X ' X ) w * tn ( X ' X ) * These two equivalence classes cannot 
be connected to each other, however, since that would require 
winding \[20 p by 2tt on only one half of the system, leading 
to an unwanted monopole in the V20 p field. This monopole 



can be removed by a phase slip, as we discuss in Section VI 



Note that there are no bulk operators local in the fermion 
variables which can distinguish the two nearly degenerate 
states. This is because all local terms on, say, wire 1 must be 
periodic in 2#i ; a term that distinguishes the two ground states 
must necessarily be odd under 0\ — > 0\ + ir. Also, our anal- 
ysis was for 2 wires. A similar analysis for N wires in series, 
coupled to the same superconducting wire, produces a degen- 
eracy of 2 A/_1 . The ground states correspond to semiclassical 
vacua +7ii7r, + n 2 ir, . . . , ^-\-tin7t), where ni = 0, 1, 
subject to the condition that the state (ni, n 2 , . . . , tin) = 
(ni + 1, n 2 + 1, . . . , n N + 1). 

As in the proximity-induced case, we can determine the 
splitting for two wires via an instanton analysis in a clean sys- 
tem with only smoothly- varying parameters, which precludes 
the possibility of phase slips. According to ( |39| ), at the soluble 
point 0+ decouples, leaving a pinned 0- field on each wire. 
To tunnel between the two vacua, we need to tunnel from 
0- = ^ to 0- = along a single length t wire, say wire 1. 
The instanton analysis proceeds exactly as that of the proxim- 
ity induced case, leading to a splitting 8E ~ exp ^— y /r K p ~ |^ . 

In such a process, 20 1 winds by ±2tt relative to V20 p . This 
can be interpreted as a vortex tunneling between the SM and 
SC wires, as depicted in process (a) in Figure]?] It can equiva- 
lently can be interpreted as a fermion tunneling from one end 
of wire 1 to the other. 

As we discuss in the next section, there is a finite region 
in parameter space around the soluble point where a similar 
analysis applies, resulting again in a splitting exponentially 
small in L 



V. PERTURBATIONS ABOUT THE SOLUBLE POINT 

We now show that the Majorana degeneracy is sta- 
ble against all possible translationally-invariant perturbations 
around the soluble point. We will consider the effect of im- 
purities and phase slips associated with them in the next sec- 
tion. The perturbations come in two varieties. First, there are 
exponentially small pair hopping terms which involve elec- 
trons in different semiconductor nano wires. Secondly, there 
are couplings between the semiconducting and superconduct- 
ing wires which we have not included in our initial model ( [37] ). 



Thirdly, there are shifts of the parameters which take us away 
from the point K\ = K 2 = 2K p , v\ = v 2 = v p . As we will 
see, the second can be accounted for with the third. 

First, we derive an effective action for interwire pair hop- 
ping, starting with the following microscopic model: 



-L/2+Z 



L/2 
L/2 



dx 1 (^\ a (x 1 )rj a (x 1 ) + h.c.) 



dx 2 (^ ] 2cT {x 2 )r] a {x 2 ) + h.c). (40) 

L/2-l 

At second order of perturbation theory in t, one obtains cross 
terms proportional to t\t 2 . These exponentially small terms 
were neglected in Eq. ( [37] ). Now we take them into account 
and study their effect on the degeneracy splitting. Consider 
the term in the Euclidean effective action proportional to t\t 2 : 



5 (12) 



2t\t 2 \ dri / dr 2 



-L/2+t 



dx\ 



L/2 



L/2 



L/2- 



(^(1)^(2)^(1)^(2) + c.c.) 



dx 2 
(41) 



(\ 2 \ 

Bosonizing the action S\ ' and integrating out the massive 
spin fields, one arrives at 



(12). 



-2tit 2 ap F 



xcos| MV-M2) |sm 



r-L/2+£ nL/2 

dr\ I dr 2 / dx\ \ dx 2 

J -L/2 JL/2-l 

'„(l)+0p(2) 



V2 



V2 



exp 



^ a /(ti-t 2 ) 2 - 



Ol^E 2 ) 2 

V 2 
F 



(2^) 2 v /(^i-^) 2 +4(n-r 2 ) 2 



■0(1)-O(2) 



(42) 



The dominant contribution to the integral over n — r 2 comes 
from short times (| T1—T2 1 <C (L—2£)/vp) and can be approx- 
imately carried out. The remaining spatial integral is peaked 
at X1/2 = ±(L/2 — £), and the action can approximately writ- 
ten as 

(12) -hh -g&(L-2l) CIPF 



S' 1Z > oc 



<pp(-Xo,T)-<j> p (x ,T) 

V2 

x sin| U-^,r)+0 p ( Xo ,r) _ g{ _^ T) _ 0{x ^ T) 



x / dr cos 



V2 



(43) 



with xq = L/2 — t. In fermionic language, the above expres- 
sion has a very clear physical interpretation: it corresponds to 
the Josephson coupling between the ends of the two wires. 
It is due to single fermion tunneling from one wire to the 
other, as depicted in process (b) in Figure [4] The terms on 
the second and third lines of ( [43] ) are bounded, so such a term 
causes a splitting which decays exponentially in (L — 2£): 

_Eg_ 

AE ~ e v f 



(L-2£) 
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We now consider interactions between the semiconductor 
and superconductor wires. We assume that because of the 
Fermi momenta mismatch in these two systems, one can ne- 
glect interactions between the charge and spin-densities at 
2k p C in the superconductor and the corresponding densi- 
ties at 2k p M in the semiconducting nanowire since these in- 
teractions will be oscillatory. We will now write down all 
possible operators couplings between the superconductor and 
the semiconductor and generate all allowed terms preserving 
U(l) symmetry. For the superconductor, the charge and spin- 
density operators are given by 



o D 



V2 

7T 

y/2. 



d x (p p 



oi = ^{Vt - ^ = - 

0*=^+^ = — cos(V26 a ) cos(V2<j> a ), 



sin 



(V26 a ) cos(V2<M 



(44) 

(45) 
(46) 
(47) 



and the singlet and triplet superconducting pairing operators 
read 



o ss = 4>\4>i - i>m 



/it,/,t 



y TS 



(l y 



_L e -iV26„ cos (^/20 o 

— e - iV * e » cos(v / 26» <T 



7ra 



-1 



sin 



(V20 a 



O. 



TS 



(48) 

(49) 

(50) 
(51) 



model p3\ one can have additional couplings which represent 
various density-density interactions 



Hi = V, 



pp 



J dxd x (j)pd x (j), 



H2 = Vg 
H3 = Vj£ 



dx 
a 

dx 



sin(v / 2^ cr ) cos(V / 2^a)^0, 
sm(V20 a ) cos(V2(i) a )d x 0. 



(58) 
(59) 
(60) 



The first term above describes the charge density-density in- 
teraction between the wires whereas the Hamiltonians in the 
second and third lines correspond to spin- spin interactions. 
The couplings between current fluctuations are similar in form 
to the density-density interactions and have been have not 
been included explicitly because their analysis is so similar. 
Assuming that \ Vaa )V ^ | are small compared to \ U\ in Eq. ( |29| ), 
the terms ( |59| ) and ( [60] ) can be dropped because the field <\> G 
orders and a is disordered. Thus, the only coupling that is 
relevant in the present setup is Hi ( [58] ). We show below that 
this quadratic term does not affect the stability of the Majo- 
rana modes. 

Therefore, a general perturbation is described by the fol- 
lowing Euclidean action: 



where SS and TS denote triplet and singlet pairing. We now 
write down these operators for the semiconductor nanowire. 
Because of the large Zeeman gap, we perform projection to 
the lowest subband as explained in Sec. [II] The charge- and 
spin-density operators in the semiconductor now become 



01 
01 

o x = 







VK 2 



a 2 p 2 F 



-dr. 



(n R -n L )-- 



d x 6 



tt ^/V 2 +a 2 p F 



-V, 



a 2 p 2 F 



(n R +n L ) = 



d x 4> 



TV 



2 ^a 2 p 2 F 



and the superconducting pairing operators read 



O 



ss 



lapp 



-a 2 p% 



a 



TS 



_ — O^s — 0. 



(52) 
(53) 
(54) 

, (55) 



(56) 
(57) 



4wles = / L " dxj dr\J^- [{d^f + viidjrf] 



7+1 



Api ; 
2tt£ 



dT[^[(d T e 2 ) 2 + vi(d x e 2 ) 2 ] 



- ^ sm(V20 p - 26 2 ) + v${d T e p ){d T e 2 ) 

J 2 l dx J dr[^ [(H) 2 + ^(Wp) 2 ] 



(61) 



The triplet pairing operators vanish because, in our model, 
the superconducting wire has a spin gap and, therefore, <j) p 
is fixed. Given these operators, one can construct all possible 
coupling terms between the superconductor and the semicon- 
ductor. In addition to the pair-hopping term, which is essen- 
tial for our proposal to work and was already included in our 



We re- write this action in terms of the new fields 9+, 0_ de- 
fined in ( 38 ). Up to local terms proportional to S(x ± ( ^ — tj ) , 
which wedrop because they will contribute negligibly to the 
bulk instanton action, we obtain 
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2 wires 



dx I dr 



A^{d T 6 + ) 2 + A ( i\d x 6 + ) 



BW(d T 0-) 2 + A£\d x 0-) 2 



2tt£ 



sin(26L) 



+ c?\d T e + ){d T e_) + c<p{d x e + ){d x oj) 



Sl, ix S 



dr 



1 ->■ 2 



2irv n 



2tt 



(62) 



where A W = Kl. + + pi, A W = 



d(i) _ A_ 

— 8ttvi 



^(1) __ 

r» W — 



K 



K p v p 



Kiv_ 



x , and similarly with 1 replaced by 2. 



27T 47T 

We see that in the general case 6+ does not decouple. How- 
ever, its action is still quadratic, so we can integrate it out ex- 
actly. We generate the following terms. On wire 1 we have 



5S^ 



duo dk 



A^^ + A^k 2 
+ J ckodk 



2(d 1} ) 2 /c 4 Ql (63) 



A^uj 2 + A x L) k 2 

We obtain an analogous expression for SS^ 2 \ We also obtain 
the following bilinear which couples wires 1 and 2: 

£S (12) = J dn dr 2 dx x dx 2 

c^cw (3 T1 a T2 (Mi)M2)}) (d T e_(i))(d T e_(2)) 
+ cpcw {d xl d X2 (e + (i)e + (2))) {a x e-(i))(d x e-(2)) 

(64) 

Here (0 + (l)0 + (2)) is the 6+ two point function between 
wires 1 and 2. 

Suppose now that the coupling SS^ were absent. Then, 
as in the instanton analysis of the proximity induced case, we 
could conclude that the lowest action instanton is translation- 
ally invariant on 1 and 2 separately; taking k — in ( [63] ) just 
gives a renormalization of the kinetic term of the 0vcenter 
of mass mode, and similarly for SS^ 2 \ We are interested in 
an instanton that tunnels from 0_ = ^ to 0_ = on wire 
1, and remains in the same vacuum on wire 2. In the present 
context, with absent, such an instanton has the same 

form as that obtained in the proximity induced case on wire 1, 
and is simply constant on wire 2. According to that analysis, 
it leads to a splitting SE ex exp(— c|). 



Now put back SS^ 12 \ We will show that the change in the 
instanton (and the change in its action) is of order j^I ' anc * 
thus negligibly small when L — 2£ ^> £. To de-clutter the 
following argument, we set all velocities equal to 1, set all 
dimensionless constants equal to 1, and let r = (x, r). The 
action is then 



S = / drdx 

J l and 2 



2 1 ^sin(20_) 



I 



d?\d?2 



+ g(r 2 - n)(d x 6-(n))(d T 6-(r 2 )) 



(65) 



where 



J{?2-n) = d Tl d T2 (fl + (fi)fl + (f 2 )> 
g(r 2 -n) = d Xl d X2 {e + {f 1 )e + (r 2 )) 



(66) 



We do not need the precise forms of / and g; rather, all we use 
is the fact that \V f(r 2 — f\)\ < ( L _^ 2 i) 3 f° r some constant 
cf, whenever x 2 — x\ > L — 2£, a condition that is always 
satisfied in ( [65] ), and a similar condition for g. Let us start 
with the instanton solution discussed above, i.e. the one that 
minimizes the action with SS^ 12 " 1 absent and tunnels between 
the two vacua only on wire 1, while staying constant in one 
of the vacua on wire 2 . We plug it into ( 65 ) and vary with 
respect to #_ (r 2 ) to obtain 



V 2 #_(f 2 ) = -^cos(2#_(f 2 )) 



ss^ 



h(r 2 ) (67) 



where 



h(f 2 ) = 7T / dfx 5(ri)d T f(f2 - ri) 



(68) 



is sourced by the instanton on 1, which, because it varies on a 
time scale £ _1 can be taken to be ttS(t\) for the purposes of 
this calculation. From (68) and the previous bound on |V/| 
we see that \h(f 2 )\ < jj^m • The key point now is that 
the dimensionful quantity h(f 2 ) is smaller than p by a factor 

of e<2 = (L-2i) 3 ^ 1 - Thus the mcms i° n of h{?2) in joTj ) 
causes 0-(f 2 ) to deviate from its zeroth order solution only 
by an amount order e. This is the first step in a perturbative 
expansion in e which shows that the inclusion of SS^ causes 
only a small change, of order e, in the instanton and its action. 

Our analysis did not require 2K p — K to be small since 
we were able to integrate out 0+ exactly regardless of their 
values. Therefore, so long as \K~ X + K~ x < 2, which im- 
plies that Ap is relevant and generates a coherence length £, 
the instanton argument is still valid and leads to a splitting 
SE oc exp(— c|). Thus, the Majorana degeneracy is stable 
over this entire region of the phase diagram, which includes 
more physically-interesting values than the soluble point. 



11 



VI. ELECTRON BACKS CATTERING AND PHASE SLIPS 

We now study the effect of processes in the superconduct- 
ing wire which backscatter a right-moving electron into a left- 
moving one or vice versa. We can include the effect of an 
electrostatic potential in the superconducting wire by adding 
a term to the action: 



pot 



dx V{x) tp1(x)ilj cr (x) 



= Jdx V(x) ^ R(J {x)^ Ra {x) + ^ L(j (x)^ L(T (x) 
+ e-^4(^ Lff (i) + e 2ik ^ La (x)^ R(T (x) 
dxV(x) \^§-d x (j>p + 2 cos(>/20p + 2k F x) cosV^o 



dxV(x) ^fd x (j) p + 2cos(V2(j) p + 2k F x) 



(69) 



In going from the penultimate line to the final one, we have 
used the fact that there is a spin gap in the SC wire which pins 
the value of <p a . The first term in the final line is harmless 
and can be absorbed by shifting (j) p which corresponds to a 
shift of the chemical potential. Therefore, we will ignore this 
term from now on. The second term in the final line causes 
2tt phase slips in the order parameter in the superconducting 



wire, e z ^ 9p , since 



y/2<t> p {x), d x (V20 p (x f ) 



— —27ri5(x — x') 



(70) 



This equation expresses the fact that when an electron in a ID 
system is backscattered, a 2tt phase slip occurs. 

These phase slips cause transitions between the two states 
of the qubit (or, in the fermion parity basis, they cause a split- 
ting between the two states). At a technical level, this occurs 
because a phase slip at the origin causes y/2Q p to wind by 2n 
on half of the system. Then 6\ can wind by 7r (while remain- 
ing at the minimum of the cosine potential), and the system 
will make a transition from the state ^f) (j, j) 

to ?f) (x> t 1 )- At a more P h y sical level > when a 
phase slip occurs, a vortex tunnels across the wire quantum- 
mechanically. Since there is no barrier for a vortex to move 
outside the wire, a vortex which tunnels through the midpoint 
of the SC wire can then encircle half of the SC wire, along 
with the NW which is in contact with that half of the SC wire. 
The vortex thereby measures the fermion parity of that NW 
by the Aharonov-Casher effect, as is depicted schematically 
in process (c) in Figure]?] 

Note that if the phase slip occurs between — L + i and L — t 
where the semiconducting wire is non-topological or is absent, 
there will only be gradient energy in 6 P (or its dual equivalent, 
fluctuation energy in (j) p ). However, if the phase slip occurs 
at a point x satisfying \L — l\ < \x\ < \L\, where there is 
a topological region of wire, then it will put a kink in y/20 p 
in a region where it is locked by the potential sin(v/2# p — 
29 1). Due to the energy cost of a kink, this will leave the 
system in a higher energy state. The kink is simply a fermion 



excited above the gap. In order to return to a ground state, 
another instanton or an anti-instanton must occur. However, 
this double process does not mix or split ground states. 

We will consider three different types of potentials V(x) 
which can backscatter electrons. First, we consider a sin- 
gle impurity. For simplicity, we will focus on the case of 
a J-function impurity at the origin, V(x) = | S(x), but the 
physics will be the same for any potential which is non-zero 
only in a region of length much less than L — 2t near the mid- 
dle of the SC wire. Then the Hamiltonian (f69b takes the form: 



imp. 



V COS 



(V20p(O)) 



(71) 



The RG equation for v follows from the scaling dimension for 

cos(V2cj) p ): 



dv 
~dl 



= (l-\K p )v 



(72) 



For K p > 2, this is irrelevant; in the large-L, low-temperature 
limit the superconductor heals itself and the backscattering 
amplitude goes asymptotically to zero. However, for K p < 2, 
the SC wire is effectively broken in two by the impurity. The 
qubit is then lost. Therefore, it is necessary to have sufficiently 
strong attractive interactions in the SC wire that K p > 2. Even 
when this is satisfied, the backscattering amplitude vanishes 
as a power-law in the system size, not exponentially. Since 
backsattering/phase slip processes cause transitions between 
the two different states of the qubit in the phase basis, they 
causes an energy splitting between states of different fermion 
parity: 



AE oc (vcos(V2</) p )) 



oc 



L k p/ 2 



(73) 



Since <p p is fixed at the ends of the SC wire (since no current 
flows off the ends), the one-point function for cos \/2(f) p has 
the L dependence show above. 

Now suppose, instead, that there is a random distribution of 
impurities so that 



V[x)V{x f ) = WS(x - x f ) 



(74) 



Then, we replicate the action by introducing an additional in- 



dex a on the field with a 



, N. We will take 



„ p - 1,2, 

N — ^ at the end of the calculation in order to take the 
quenched average over all realizations of the disorder. The 
disorder-averaged effective action takes the form: 

Random = j dr dr' dx W cos ( V2 [<t>« (x, r) - rf(x,r f )]) 

(75) 

The RG equation for W is: 

dW 



dl 



= (3 - K p ) W 



(76) 



Thus, we need a larger K p for the superconductivity to survive 
a random distribution of impurities, and if K p > 3 is satisfied, 
then there will be an energy splitting: 



AE oc 



W 



L K p~ 



(77) 
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FIG. 4: A schematic depiction of the different processes which split 
the two states of the qubit. (a) A vortex can tunnel between a semi- 
conducting wire and the superconducting wire. This causes a split- 
ting which is exponentially- small in £. (b) An electron can tunnel 
from one semiconducting wire to another through the superconduct- 
ing wire. This causes a splitting which is exponentially- small in 
L — 2L (c) A vortex can tunnel through the superconducting wire be- 
tween the semiconducting wires as a result of an electron backscat- 
tering process. This causes a splitting which decays as a power of 
L. This model applies also to a situation in which there is a single 
semiconducting wire which extends from —L/2 to L/2 and is in a 
topological phase for L/2 — £ < \x\ < L (similar to our two SM 
wires) and is in a non-topological phase for \x\ < L/2 — i (similar 
to our superconducting wire). 



Thus far, we have focused on backscattering by impurities, 
which effectively create weak spots in the wire where a vortex 
can tunnel through. However, even in a completely clean sys- 
tem, there is some amplitude for backscattering. For instance, 
let us suppose that V(x) is constant near the middle of the 
wire and goes to zero smoothly near the ends. To make this 
concrete, let us take V(x) — Vo for \x\ «L/2 and V(x) =0 
for \x\ = L/2. We will assume that V(x) varies smoothly, 
so that the Fourier transform V(q) ~ e~ q b forg^>l/L, 



where b > £. Then, from the last line of Eq. 69 we expect a 
splitting 



AE oc 



/ 



dx 



V(x) cos2kpx 



< 



-+£ 



(f-l; 



(78) 



Therefore, as the potential becomes smoother and smoother, 
the splitting which it induces through electron backscatter- 
ing/phase slips goes exponentially to zero with the length 
scale b over which the potential varies. Inhomogeneities en- 



hance backscattering, as we saw in Eqs. 73 77 



VII. DISCUSSION 



As we saw in the Section [VlJ the effects of electron 
backscattering by impurities can be mitigated by making K p 
large. In a superconducting wire, K p 

k 2 F A w CX Channels and V 

being the cross-sectional area, superconducting stiffness and 
compressibility, respectively 28 . Therefore, if the supercon- 
ducting wire has enough channels or, equivalently, if the su- 
perconducting wire has sufficiently large cross-sectional area 



27TV 'A w p s K CX 

\J A w p s /K with A w , p s and n 



and/or sufficiently large superfluid density, we can have K p 
large. Although this is not as good as exponential decay as a 
matter of principle, it may be just as good as a practical mat- 
ter. This may be important since it could be very difficult to 
tune the chemical potential appreciably in the semiconducting 
wire (which is necessary to move the Majorana zero modes) 
if it is in contact with a bulk 3D superconductor. Furthermore, 
coating the semiconducting wires with superconducting ma- 
terial, as depicted in Figure [2] may be the easiest way to make 
a complex network of wires (especially a three-dimensional 
network) which is in contact with a superconductor 47 . How- 
ever, such an architecture will necessarily be, at best, an 
algebraically-ordered superconductor (except, perhaps, at the 
lowest temperatures, at which the coupling between wires 
causes a crossover to 3D superconductivity). Therefore, it is 
significant that our results show that such a network supports 
Majorana fermion nearly-zero modes and that their splitting 
can be made small (albeit not exponentially so). 

We also note that it is only important that K p be large in 
the regions between the topological semiconducting wire seg- 
ments. In the topological semiconducting wire segments, the 
phase is locked so that 2tt phase slips cannot occur (although 
harmless 47r phase slips can occur). Therefore, one can imag- 
ine a scenario in which the topological segments are coated 
with a thin superconducting film while the non-topological 
segments between them are in contact with essentially bulk 
3D superconductors. This would lead to a protected topo- 
logical qubit, although it would be difficult (if not impossi- 
ble) to move the Majorana zero modes since that would in- 
volve tuning the chemical potential in the non-topological re- 
gions - which are in contact with bulk 3D superconductors 
- to drive them into the topological phase. One may al- 
ternatively, in a system in which a semiconducting wire is 
coated with a thin layer of superconducting material, use a 
gate voltage to occupy a large (even) number of sub-bands 
of the semiconducting wire in the non-topological regions. 
This would lead to a large effective K p for the combined 
superconductor- semiconductor system in the non-topological 
regions and, therefore, a large power for the decay of the split- 
ting due to phase slips in these regions. 

As the previous sentence anticipates, our methods should 
be generalizable to multi-channel semiconducting wirePSEIl 
They should also apply to a semiconducting wire which is 
near a superconducting grain (as in a model of quasi- ID wires 
in LAO/STO interfaces 24 ). If the linear size of the grain, r, 
is smaller than the superconducting coherence length, £, then 
we can treat the grain as a zero-dimensional system. Suppose 
that the wire also has length r. Then the Hamiltonian for the 
wire coupled to the grain is simply ([T]) with cj> independent of 
position i but dependent on time. There will also be a charg- 
ing energy U(N — Nq) 2 which causes <p to fluctuate. There 
will be no long-ranged order in the superconducting grain, but 
it can still induce a single-fermion gap in the semiconduct- 
ing wire. Of course, if the wire has length L ^> £, then the 
grain will only change the behavior of a short section of the 
wire, and the two ends of this section will be relatively close to 
each other. But if the wire passes near many such grains, then 
they can induce a single-fermion gap in the wire. If the cou- 
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pling between the grains is large compared to their charging 
energies, then, in the long- wavelength limit, the grains will 
develop algebraic order. The superconducting grains can be 
modeled by a superconducting wire, and this situation can be 



modeled with the Hamiltonian of Section III but with a very 
small velocit}JSI. If there is ohmic dissipation, then the grains 
may not even have power-law superconducting order but may 
have exponentially-decaying superconducting correlations. 

In fact, we will have Majorana zero modes in a system 
with exponentially-decaying superconducting correlations if 
we simply take our model to finite temperature. Then, the 
6+ field in Eq. ([39]) will have exponentially-decaying corre- 
lations, with a correlation length inversely proportional to the 
temperature. The 6- field will still be pinned to a minimum 
of the potential, but it will be possible for the system to be 
thermally-excited over the barrier from one minimum to the 
other. Therefore, if is the bulk single-fermion gap, there 
will be a contribution due to processes (a) and (b) in Fig. |4]to 
the coherence time for a Majorana qubit of order ~ e~ A ^ T , 
just as if there were long-ranged superconducting order. How- 
ever, there will also be a contribution from quantum phase 
slips, process (c), which will increase with temperature as 
T Kp l 2 for a single impurity and T Kp ~ 2 for a random distri- 
bution of impurities. We similarly expect Majorana fermion 
zero modes to survive in two-dimensional structures in which 
a superconducting gap is induced via the proximity effect to 
stabilize a phase with Ising anyond^Elbut long-ranged super- 
conducting order is disordered by quantum or thermal fluctu- 
ations. If the single-particle gap remains, then the Majorana 



fermion zero modes associated with the Ising anyons could 
survive. However, quantum phase slips are suppressed 40 ^ and, 
therefore, the splitting will be exponentially, rather than al- 
gebraically, decaying. Of course, there is nothing surprising 
about having protected Majorana zero modes in a system with 
no long-range order or even algebraic order, since this is pre- 
cisely the case with any true topological phase of matter, as 
in the examples mentioned in the introduction. However, the 
particular route which we have found to such a system is new 
and interesting. 

In this paper, we have shown that a gapless system can be 
nearly as good as a fully-gapped one at supporting protected 
Majorana fermion zero modes. It is an interesting open ques- 
tion whether a gapless system might be capable of support- 
ing protected degrees of freedom which cannot occur in fully- 
gapped ID systemsS^ED. 

Note added: After the initial version of this paper appeared 
on the arXiv, several other paperJ^S^ on related topics were 
submitted to the arXiv. 
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